Superconducting single-charge transistor in a tunable dissipative environment 
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We study a superconducting single-charge transistor, where the coherence of Cooper pair tunneling 
is destroyed by the coupling to a tunable dissipative environment. Sequential tunneling and cotun- 
neling processes are analyzed to construct the shape of the conductance peaks and their dependence 
on the dissipation and temperature. Unexpected features are found due to a cross-over between 
two distinct regimes, one 'environment-assisted' the other 'environment-dominated'. Several of the 
predictions have been confirmed by recent experiments. The model and results apply also to the 
dynamics of Josephson junction quantum bits on a conducting ground plane, thus explaining the 
influence of dissipation on the coherence. 
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The behavior of a quantum system coupled to a dis- 
sipative environment is one of the paradigms of modern 
physics |l|]. It is the central challenge in any attempt 
to build and manipulate quantum information systems, 
since fluctuations and dissipation limit the quantum co- 
herence of the device. For instance, for qubits based on 
two different charge states of a Josephson junction g] 
one important source of dissipation are the normal elec- 
trons in the ground plane. A conducting ground plane is 
an appropriate tool however, to compensate for the ran- 
dom offset charges of the superconducting island. In the 
experimental realization of Nakamura et al. H this was 
achieved by mounting the device on a gold ground plane. 

In earlier experiments the Berkeley group B had 
mounted an array of Josephson junctions on a two- 
dimensional electron gas (2DEG), separated from it by 
an insulating layer. In this setup the density of nor- 
mal electrons in the 2DEG, and hence the source and 
strength of the dissipation, can be tuned by a back gate 
over a wide range (see also S). This influences the col- 
lective properties of the array in an qualitative way [|6| . 
In more recent experiments [[7| the Berkeley group in- 
vestigated the transport properties of a Superconducting 
Single-Charge Transistor (SSCT) coupled to a tunable 
2DEG (see Fig. [^). They observed an unexpected de- 
pendence of the conductance peak heights and widths on 
temperature and dissipation. 

In previous work the nonlinear current branch of a 
SSCT has been studied |p-[l0||. Due to the couphng of 
the SSCT to the metallic ground plane the Cooper pair 
charges interact with normal metal image charges, which 
turns the tunneling of Cooper pairs dissipative. A cen- 
tral result of the present paper is the identification and 
characterization of the previously uninvestigated nor- 
mal current branch in the linear and nonlinear response 
regime. We determine how the conductance depends on 
the gate and transport voltage, environment conductance 
and temperature. We also provide estimates how dissi- 
pation limits the quantum oscillations of qubits. 
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FIG. 1. (a) Top view of the system under consideration. 
The 2DEG is in the substrate, (b) Equivalent circuit. 



In the following we concentrate on systems (see Fig. |l|) 
with low resistance per square of the ground plane, R <C 
/i/4e^ « 6.5 kri, and we choose the gate and junction 
capacitances, Cg and Cj, and the capacitances between 
the island and the leads to the metallic plane, Cq and C\, 
to be ordered as follows: Cg, Cj <C Co <C Ci. The effects 
of the dissipative metallic ground plane on the SSCT are 
captured by the real part of the impedance. 



Re{Z{uj)} = R/[l + (LuRCjf] . 



(1) 



Interestingly, even though the different capacitances in- 
troduce several frequency regimes, in all regimes TieZ{(jj) 
is determined by the capacitance of the tunnel junction 
Cj only. Quasiparticle dissipation may be ignored at tem- 
peratures much below the superconducting gap. 

We construct our theory following Ref. |jri|,|l2|. The 
environment can be described by a harmonic oscillator 
bath. An integration over the quadratic bath degrees of 
freedom is possible. The resulting sequential tunneling 
rate through one junction |l3|. 



T{SE,i,) = {^/2h)EJP{6E,^) 



(2) 



depends on the change in charging energy during the tun- 
neling process, which in turn depends on the gate voltage 
ng = VgCg/e and the transport voltage Vtv For a process 
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FIG. 2. Coiraactance peaks in lowest order (^iquential tun- 
neling) scaled by the critical current of a single junction, Ic 
and the charging energy Ec of the island, assuming Cq <^ Cj . 



which increases the number of excess Cooper pairs from 
71 to n -|- 1 is given by 



SE,Y, = -AEci2n - uq + 1) + eVt, 



(3) 



The (total) capacitance of the island determines the 
energy scale Eq = e^/2(2Cj + Cg), while Ej is 
the Josephson coupling energy of the tunnel junc- 
tions (here assumed to be equal). The function P{E) 
can be expressed by a correlation function K{t) = 
¥ I-oo ^^^{Z{^)} {coth (1^) [cos(^i) - 1] - isin(c.O} 



In the situation considered, K(t) can be evaluated an- 
alytically and in the long-time limit, |t| ^ -^Cj, it reads 



i^(t) = --Lr|f|+logfl-e-2'^^l*l 
L V 



7 + log(27r™Cj)-Hzsign(t)- 



In the low energy regime, i?, T <C RCj, wc find 
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(5) 



where 7 = 0.577..., g = h/Ae^R is the dimensionlcss 
conductance of the ground plane, and B{x, y) the Beta 
function. We observe that the energy appears only in the 
dimensionlcss combination e — Eg /2'kT . 

Given the tunneling rates T{5Ec\i) we construct a mas- 
ter equation and determine the current / through the 
SSCT. The difference of the charging energies SEch varies 
as the gate voltage is tuned. At low temperatures, for 



most values of rig, the large cost in charging energy sup- 
presses transport by the Coulomb blockade. At the de- 
generacy points of the states with n and n -t- 1 Cooper 
pairs (i.e. dEch = 0), transport becomes possible. Thus 
the conductance exhibits Coulomb oscillations as a func- 
tion of the gate voltage. 

Previously published work |13,M concentrated on the 
low temperature regime e 3> 1. Here the asymptotic ex- 
pansion of Eq. (||) yields P{E) « 2/gE [0. For instance 
at the degeneracy point, we have 5Ec\i — eVtr and the I- 
V characteristics is strongly nonlinear 



/ = AnE^^/hgVt, . 



(6) 



Away from the degeneracy point the high energy cost, 
SEch ^ T, suppresses tunneling (Coulomb blockade). 
Thus this regime is still described by |e| ^ 1. Detailed 
balance requires P{—E) = e~^/'^P{E), from which one 
finds that the conductance is suppressed by an exponen- 
tial activation factor ||llf| . 

A closer inspection shows, however, that for typical 
transport voltages of Vtr ~ 1 nV, even at the lowest 
temperatures of T = 20mK, for the typical values of 
g = 1..100 one has e <C 1. This motivated us to study this 
previously ignored parameter regime. Surprisingly we 
find that both at the peaks and in the Coulomb blockade 
regime the system's response becomes linear, although 
for different reasons. 

First, at the peaks and for low transport voltages 
SEch = eVtr <C T one can expand Eq. (||). We find that 
in this regime the current response is linear 



/(Ftr, riG = n+l) = eTTE]geVt,/hT^ 



(7) 



and the height of the conductance peak now increases 
with increasing conductance g of the metallic plane, in 
contrast to what could be concluded from Eq. (^. 

The full shape of the conductance peak derived from 
Eq. (^ is shown in Figure 2. At low enough transport 
voltages, Vtr < 27rT/g, at the peaks we have e < 1 and 
the current response is linear and transport is provided 
by ' environment- rfommafed' sequential tunneling (ED- 
ST). The peak height is proportional to 1/T^ at fixed 
g (Fig. 11(a)) and proportional to g at fixed T (Fig. ||(b)). 
Away from the peak the energy change 6Ech for a tun- 
neling process is large even for Vtr = 0. Therefore one 
crosses over to the regime e 3> 1. Here the differen- 
tial conductance at some fixed voltage is proportional to 
l/g. This tendency is opposite to that of the peaks, and 
manifests itself in the crossing of the conductance G{Vg) 
curves in the tails (see Fig. 0(b)). In this 'environment- 
assisted' sequential tunneling (EA-ST) regime, the con- 
ductance may still be appreciable; it is reduced due to 
Coulomb blockade only for SE^h > T, i.e. at e > g. 
We also find that the width of the peak increases with 
decreasing g, and the conductance in the tails decreases 
with increasing T. 



In a transport experiment the voltage is necessarily fi- 
nite, Vtr > 0. At low enough temperatures T <C Vtrff/STr 
this forces the system to e > 1, where the nonlinear re- 
sponse formulae govern the physics. Thus, the apparent 
divergence of the conductance at the peak, G ~ 1/T^ 
is regularized by the transition to the nonlinear regime, 
yielding a finite conductance for T — > (Fig. ^(c)). Also, 
at sufficiently large Vtr the conductance monotonically 
increases with decreasing g over the whole energy range 
(Fig. 1(d)). 

In second order in Ej/Ech the current is exponentially 
suppressed in the Coulomb blockade regime. Here, higher 
order terms, such as the fourth order cotunneling pro- 
cess, may yield important contributions. Cotunneling de- 
scribes the simultaneous coherent transfer of two Cooper 
pairs, where the the energetically forbidden intermediate 
state is occupied only virtually. As a consequence the 
process is not suppressed by an exponential activation 
factor in the blockade regime. The most relevant term 
(for SEch < -T/g) leads to the rate [^ 
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At low voltages we thus find again a linear response 
regime with a decay rate Fct oc Vtrg/T'^- 

Combining these results, we can read off the peak 
width within linear response. Depending on the ex- 
perimental conditions, the conductance is dominated 
by environment-assisted sequential tunneling (EA-ST) 
or environment-dominated cotunneling (ED-CT). This 
leads to a half-width of i5£'ch,i/2 oc max(r/g,_Ej) 
whereas a fixed off-peak conductance go is reached at 

SEch,go oc max f^j/^/ggo, \/g/go Ej/Tj. At high trans- 
port voltages and large SEch the response becomes non- 
linear, yielding Fct oc 1/gV^. Remarkably however the 
ED-ST regime extents to large Vtr, in the region where 
the bare charging energy difference and Vtr combine to a 
small net SEch, generating again a low value of |e| <C 1. 
Here again the incoherent transport of Cooper pairs dom- 
inates, with a conductance G oc gE^j/T^. The expres- 
sions for the conductances in the different parameter 
regimes are summarized in Fig. 3. 

So far we considered relatively short and wide leads, 
when the bilayer formed by the leads and the metallic 
plane can be described as a capacitor. If the leads are 
longer and narrower, the spatial correlations of charges 
moving within the bilayer modify the fluctuation spec- 
truma and should be taken into account. Then the proper 
model for the leads is an RC transmission line using the 
per-square capacitance and resistance, rather than the 
CL-capacitors in Fig. 1. The impedance of such lines is 
Zrc{^) = \/R/ii^Gsq where Cgq is the square capaci- 
tance between the lead and the ground plane. In this 
case the aj~^/^ divergence governs the impedance. At 
T = the model has been analysed before M,f2| . Here 




FIG. 3. Overview of the different transport regimes as a 
function of the transport voltage Vtr and the gate voltage (en- 
tering in the combination 5Ech,o = —{n — no — l/2)4e'^/C). 
ED-CT: Environment-dominated cotunneling; NL-CT: Non- 
linear cotunneling; EA-ST: Environment-assisted sequential 
tunneling; ED-ST: Environment-dominated sequential tun- 
neling; NL-ST: Nonlinear sequential tunneling 



we concentrate on the high temperature limit. At the 
peak and for T }^ E we expand in powers of UJ/2T to 
obtain 

K{t) = ~2^7r\t\/ gCsc, (4|i|T - isign(i)) 

which is dominated by the |ip/^ term. Then 
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The response of the system is again linear, but the con- 

dutance scales with new exponents: G oc T~^/^ and 
Gocgi/3. 

If the SSCT is driven into the normal state, its own 
dynamics becomes dissipative, but it is still strongly in- 
fluenced by the environment. In contrats to the resonant 
tunneling situation of the superconducting state we have 
to sum now over the normal flnal states. In this case the 
tunneling rate is found from P{E), cq. a along the lines 

of El 
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It depends on the energy again through e. Since the 
charge carriers are now normal electrons, g is redefined 
accordingly as g = h/e^R and the junction strength is 
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FIG. 4. Cofefetance peaks in lowest orde^fe^^erturbation 
theory for the normal state, scaled by the tunneling conduc- 
tance Gt of the junctions. 

characterized by the tunneling resistance Rt- For e <C 1 
we find Tl = {T/Rt){l + SE^h/T) for g -^ oo. Thus at 
the peak the response is linear again, but this time with 
a temperature independent value as T ^ 0. For e 3> 1, 
the rate displays the usual Coulomb blockade 

r(5^ch) = d{SE^^) ^^'^^ , J e-^SE^i, RCtf/a 

V ch; V '='^^7ri?tr(2-|-2^)^ ' 

The change of the power-laws can be understood by the 
additional phase space factor ex E"^ /T"^ in the normal 
state. The resulting independence of the environment 
is caused by the fact that the system is dissipative by 
itself and the external dissipation is weak at large g. Our 
numerical results in Fig. show the (weak) dependence 
on the environment conductance in the T — 0-limit. 

We are now ready to review the experimental results 
of the Berkeley group 0. They constructed the SSCT 
just like that of Fig. 1. The parameters were: Ej = 1 
K, Ct = 0.3fF,Co = 1.5fF. The backplane conductance 
was swept in the range g = 1...100 and the tempera- 
ture in the range T — 20. ..1000 mK. They found in the 
Coulomb oscillations of the conductance a linear response 
regime at the peaks, confirming a major prediction of 
the present theory. The peak conductance was found 
to scale as G ex g°'^ and G ex T~°-^. Their lead arrange- 
ment is modeled more adequately with a transmission 
line. Our relevant results for this case are: G oc g^^^ and 
G ex T~^/^. When the SSCT was driven normal by a 
large magnetic field, they found that the peak conduc- 
tance exhibits very weak dependence on g and T. Our 
theory predicted that G goes to a constant for either 
large g or small T. So several of the measured g and T 
dependences are in agreement with the present theory. 
Finally, between the peaks they do not find an exponen- 
tial suppression of G, making it plausible that indeed the 
cotunneling process is the dominant transport channel, 
as described above. 

A similar system, a superconducting electron box 
mounted on a metallic ground plane [16| , has been used to 
demonstrate the coherent oscillations of charge states [|| , 
which makes this setup a promising candidate for a qubit 
in quantum computers. To this end, it has to be oper- 
ated at the peak. The quality factor Q of the oscillations 
can be determined by standard methods M. At the low 



temperatures oi T <S^ Ej we obtain 

Q = cot(7r/(.g - 2)) w g/n 



(8) 



as for the metallic case, g ^ 2. For the parameters of S] 
we find Q in the range of several hundreds. Thus the dis- 
sipation allows of the order of a thousand cycles of the 
quantum oscillations. While this is encouraging, it still 
shows that the metal planes, which were introduced to 
screen the random charges, induce a potentially serious 
limitation on the maximum number of coherent opera- 
tions. 
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